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We report on a search for gravitational waves from binary black hole inspirals in the data from the second 
science run of the LIGO interferometers. The search focused on binary systems with component masses between 
3 and 20 Mq . Optimally oriented binaries with distances up to 1 Mpc could be detected with efficiency of at 
least 90%. We found no events that could be identified as gravitational waves in the 385.6 hours of data that we 
searched. 
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I. INTRODUCTION 

The Laser Interferometric Gravitational Wave Observatory 
(LIGO) Ql consists of three Fabry-Perot-Michelson interfer- 
ometers, which are sensitive to the minute changes that would 
be induced in the relative lengths of their orthogonal arms by a 
passing gravitational wave. These interferometers are nearing 
the end of their commissioning phase and were close to design 
sensitivity as of March 2005. During the four science runs 
that have been completed until now (first (SI) during 2002, 
second (S2) and third (S3) during 2003 and fourth (S4) dur- 
ing 2005) all three LIGO interferometers were operated stably 
and in coincidence. Although these science runs were per- 
formed during the commissioning phase they each represent 
the best broad-band sensitivity to gravitational waves that had 
been achieved up to that date. 

In this paper we report the results of a search for gravi- 
tational waves from the inspiral phase of stellar mass binary 
black hole (BBH) systems, using the data from the second 
science run of the LIGO interferometers. These BBH sys- 
tems are expected to emit gravitational waves at frequencies 
detectable by LIGO during the final stages of inspiral (decay 
of the orbit due to energy radiated as gravitational waves), 
the merger (rapid infall) and the subsequent ringdown of the 
quasi-normal modes of the resulting single black hole. 

The rate of BBH coalescences in the Universe is highly un- 
certain. In contrast to searches for gravitational waves from 
the inspiral phase of binary neutron star (BNS) systems 01, 
it is not possible to set a reliable upper limit on astrophys- 
ical BBH coalescences. That is because the distribution of 
the sources in space, in the component mass space and in the 



* Currently at Intel Corp. 

* Currently at University of Tours, France 
&& Currently at Lightconnect Inc. 
Currently at W.M. Keck Observatory 

Currently at ESA Science and Technology Center 
++ Currently at Raytheon Corporation 

^ Currently at New Mexico Institute of Mining and Technology / Magdalena 

Ridge Observatory Interferometer 

** Currently at Mission Research Corporation 

TT Currently at Harvard University 

H Currently at Lockheed-Martin Corporation 

Permanent Address: University of Tokyo, Institute for Cosmic Ray Re- 
search 

" " Permanent Address: University College Dublin 

II I' Currently at University of Alaska Anchorage 
Currently at Research Electro-Optics Inc. 

"Currently at Institute of Advanced Physics, Baton Rouge, LA 

*** Currently at Thirty Meter Telescope Project at Caltech 

+ + + Currently at European Commission, DG Research, Brussels, Belgium 

*** Currently at University of Chicago 

## Currently at LightBit Corporation 

^Permanent Address: IBM Canada Ltd. 

tt Currently at The University of Tokyo 

§§§ Currently at University of Delaware 

III Permanent Address: Jet Propulsion Laboratory 

II II II Currently at Shanghai Astronomical Observatory 
Currently at Laser Zentmm Hannover 



spin angular momentum space is not reliably known. Addi- 
tionally, the gravitational waveforms for the inspiral phase of 
stellar-mass BBH systems which merge in the frequency band 
of the LIGO interferometers are not known with precision. We 
perform a search that aims at detection of BBH inspirals. In 
the absence of a detection, we use a specific nominal model 
for the BBH population in the Universe and the gravitational 
waveforms given in the literature to calculate an upper limit 
for the rate of BBH coalescences. 

The rest of the paper is organized as follows. Sec.llTlpro- 
vides a short description of the data that was used for the 
search. In Sec.[II1]we discuss the target sources of the search 
and we explain the motivation for using a family of phe- 
nomenological templates to search the data. In Sec. I1VI we 
give a detailed discussion of the templates and the filtering 
methods. In Sec. IV A| we provide information on various data 
quality checks that we performed, in Sec. IV Bl we describe in 
detail the analysis method that we used and in Sec. IV Cl we 
provide details on the parameter tuning. In Sec. I VII we de- 
scribe the estimation of the background and in Sec. I VIII we 
present the results of the search. We finally show the calcula- 
tion of the rate upper limit on BBH coalescences in Sec. I Villi 
and we provide a brief summary of the results in Sec. lIXI 



II. DATA SAMPLE 

During the second science run, the three LIGO interferome- 
ters were operating in science mode (see Sec. IV Al l. The three 
interferometers are based at two observatories. We refer to the 
observatory at Livingston, LA, as LLO and the observatory at 
Hanford, WA as LHO A total of 536 hours of data from the 
LLO 4 km interferometer (hereafter LI), 1044 hours of data 
from the LHO 4 km (hereafter HI) interferometer, and 822 
hours of data from the LHO 2 km (hereafter H2) interferom- 
eter was obtained. The data was subjected to several quality 
checks. In this search, we used only data from times when 
the LI interferometer was running in coincidence with at least 
one of HI and H2, and we only used continuous data of dura- 
tion longer than 2048 s (see Sec. lVBk After the data quality 
cuts, there was a total of 101.7 hours of LI -HI double coinci- 
dent data (when both LI and HI but not H2 were operating), 
33.3 hours of L1-H2 double coincident data (when both LI 
and H2 but not HI were operating) and 250.6 hours of Ll- 
H1-H2 triple coincident data (when all three interferometers 
were operating) from the S2 data set, for a total of 385.6 hours 
of data. 

A fraction (approximately 9%) of this data (chosen to be 
representative of the whole run) was set aside as "playground" 
data where the various parameters of the analysis could be 
tuned and where vetoes effective in eliminating spurious noise 
events could be identified. The fact that the tuning was per- 
formed using this subset of data does not exclude the possibil- 
ity that a detection could be made in this subset. However, to 
avoid biasing the upper limit, those times were excluded from 
the upper limit calculation. 

As with earlier analyses of LIGO data, the output of the an- 
tisymmetric port of the interferometer was calibrated to obtain 
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a measure of the relative strain AL/L of the interferometer 
arms, where AL = L x — L y is the difference in length be- 
tween the x arm and the y arm and L is the average arm length. 
The calibration was measured by applying known forces to the 
end mirrors of the interferometers before, after and occasion- 
ally during the science run. In the frequency band between 
100 Hz and 1500 Hz, the calibration accuracy was within 10% 
in amplitude and 10° of phase. 



III. TARGET SOURCES 

The target sources for the search described in this paper 
are binary systems that consist of two black holes with com- 
ponent masses between 3 and 20 Mq, in the last seconds 
before coalescence. Coalescences of binary systems consist 
of three phases: the inspiral, the merger and the ringdown. 
We performed the search by matched filtering the data us- 
ing templates for the inspiral phase of the evolution of the 
binaries. The exact duration of the inspiral signal depends on 
the masses of the binary. Given the low-frequency cutoff of 
100 Hz that needed to be imposed on the data (see Sec. lVBf 
the expected duration of the inspiral signals in the S2 LIGO 
band as predicted by post-Newtonian calculations varies from 
0.607 s for a 3 - 3 M© binary to 0.013 s for a 20 - 20 M 
binary. 

The gravitational wave signal is dominated by the merger 
phase which potentially may be computed using numerical 
solutions to Einstein's equations. Searching exclusively for 
the merger using matched-filter techniques is not appropri- 
ate until the merger waveforms are known. BBH mergers are 
usually searched for by using techniques developed for de- 
tection of unmodeled gravitational wave bursts j3[. However, 
for reasons that will be explained below, it is possible that the 
search described in this paper was also sensitive to at least 
part of the merger of the BBH systems of interest. Certain 
re-summation techniques have been applied to model the late 
time evolution of BBH systems which makes it possible to 
evolve those systems beyond the inspiral and into the merger 
phase B II S IE and the templates that we used 

for matched filtering incorporate the early merger features (in 
addition to the inspiral phase) of those waveforms. 

The frequencies of the ringdown radiation from BBH sys- 
tems with component masses between 3 and 20 Mq range 
from 295 Hz to 1966 Hz (U Q2 H and the gravitational 
wave forms are known. Based on the frequencies of these sig- 
nals, some of the signals are in the S2 LIGO frequency band 
of good sensitivity and some are not. At the time of the search 
presented in this paper, the matched-filtering tools necessary 
to search for the ringdown phase of BBH were being devel- 
oped. In future searches we will look for ringdown signals 
associated with inspiral candidates. 

Finally, we have verified through simulations that the pres- 
ence of the merger and the ringdown phases of the gravita- 
tional wave signal in the data does not degrade our ability to 
detect the inspiral phase, when we use matched filter tech- 
niques. 



A. Characteristics of BNS and BBH inspirals 

We use the standard convention c = G = 1 in the remainder 
of this paper. 

The standard approach to solving the BBH evolution prob- 
lem uses the post-Newtonian (PN) expansion 1 10] of the Ein- 
stein equations to compute the binding energy E of the binary 
and the flux F of the radiation at infinity, both as series ex- 
pansions in the invariant velocity v (or the orbital frequency) 
of the system. This is supplemented with the energy balance 
equation (dE/dt = —F) which in turn gives the evolution 
of the orbital phase and hence the gravitational wave phase 
which, to the dominant order, is twice the orbital phase. This 
method works well when the velocities in the system are much 
smaller compared to the speed of light, d<1, Moreover, the 
post-Newtonian expansion is now complete to order v 7 giv- 
ing us the dynamics and orbital phasing to a high accuracy 
fl4l fl5ll . Whether the waveform predicted by the model to 
such high orders in the post-Newtonian expansion is reliable 
for use as a matched filter depends on how relativistic the sys- 
tem is in the LIGO band. For the second science run of LIGO, 
the interferometers had very good sensitivity between 100 and 
800 Hz so we calculate how relativistic BNS and BBH sys- 
tems are at those two frequencies. 

The velocity in a binary system of total mass M is related 
to the frequency / of the gravitational waves by 

v = (ttM/) 1 / 3 . (1) 

When a BNS system that consists of two 1.4 Mq compo- 
nents enters the S2 LIGO band, the velocity in the system is 
v ~ 0.16; when it leaves the S2 LIGO band at 800 Hz, it is 
v ~ 0.33 and the system is mildly relativistic. Thus, relativis- 
tic corrections are not too important for the inspiral phase of 
BNS. 

BBH systems of high mass, however, would be quite rela- 
tivistic in the S2 LIGO band. For instance, when a 10— 10 Mq 
BBH enters the S2 LIGO band the velocity would be v ~ 
0.31. At a frequency of 200 Hz (smaller than the frequency 
of the innermost stable circular orbit, explained below, which 
is 220 Hz according to the test-mass approximation) the ve- 
locity would be v ~ 0.40. Such a binary is expected to 
merge producing gravitational waves within the LIGO fre- 
quency band. Therefore LIGO would observe BBH systems 
in the most non-linear regime of their evolution and thereby 
witness highly relativistic phenomena for which the perturba- 
tive expansion is unreliable. 

Numerical relativity is not yet in a position to fully solve 
the late time phasing of BBH systems. For this reason, in 
recent years, (non-perturbative) analytical resummation tech- 
niques of the post-Newtonian series have been developed to 
speed up its convergence and obtain information on the late 
stages of the inspiral and the merger 11611 . These resummation 
techniques have been applied to the post-Newtonian expanded 
conservative and non-conservative part of the dynamics and 
are called effective-one-body (EOB) and P-approximants 
(also referred to as Pade approximants) I^SSSSIal- Some 
insights into the merger problem have been also provided in 
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lll7tll8ll by combining numerical and perturbative approxima- 
tion schemes. 

The amplitude and the phase of the standard post- 
Newtonian (TaylorT3, 0), EOB, and Pade waveforms, eval- 
uated at different post-Newtonian orders, differ from each 
other in the last stages of inspiral, close to the innermost sta- 
ble circular orbit (ISCO, fl6lo . The TaylorT3 and Pade wave- 
forms are derived assuming that the two black holes move 
along a quasi-stationary sequence of circular orbits. The EOB 
waveforms, extending beyond the ISCO, contain features of 
the merger dynamics. All those model-based waveforms are 
characterized by different ending frequencies. For the quasi- 
stationary two-body models the ending frequency is deter- 
mined by the minimum of the energy. For the models that 
extend beyond the ISCO, the ending frequency is fixed by the 
light-ring Ki ll 911 of the two-body dynamics. 

We could construct matched filters using waveforms from 
each of these families to search for BBH inspirals but yet the 
true gravitational wave signal might be "in between" the mod- 
els we search for. In order not to miss the true gravitational 
wave signal it is desirable to search a space that encompasses 
all the different families and to also search the space "in be- 
tween" them. 



B. Scope of the search 

Recent work by Buonanno, Chen and Vallisneri |l9| (here- 
after BCV) has unified the different approximation schemes 
into one family of phenomenological waveforms by introduc- 
ing two new parameters, one of which is an amplitude cor- 
rection factor and the other a variable frequency cutoff, in 
order to model the different post-Newtonian approximations 
and their variations. Additionally, in order to achieve high 
signal-matching performance, they introduced unphysical pa- 
rameters in the phase evolution of the waveform. 

In this work we used a specific implementation of the phe- 
nomenological templates. As these phenomenological wave- 
forms are not guaranteed to have a good overlap with the true 
gravitational wave signal it is less meaningful to set upper 
limits on either the strength of gravitational waves observed 
during our search or on the coalescence rate of BBH in the 
Universe than it was for the BNS search in the S2 data 01. 
However, in order to give an interpretation of the result of 
our search, we did calculate an upper limit on the coalescence 
rate of BBH systems, based on two assumptions: (1) that the 
model-based waveforms that exist in the literature have good 
overlap with a true gravitational wave signal and (2) that the 
phenomenological templates used have a good overlap with 
the majority of the model-based BBH inspiral waveforms pro- 
posed in the literature lll9ll . 

To set the stage for later discussion we plot in Fig. ^the 
distance at which a binary of two components of equal mass 
that is optimally oriented (positioned directly above the in- 
terferometer and with its orbital plane perpendicular to the 
line of sight from the interferometer to the binary) would pro- 
duce a signal-to-noise ratio (SNR, see Sec. I1VI of 8 in the 
LIGO interferometers during the second science run. We re- 



fer to this distance as "range" of the interferometers. The solid 
line shows the range of the LIGO interferometers for matched 
filtering performed with the standard post-Newtonian (Tay- 
lorT3) waveforms, which predict the evolution of the system 
up to the ISCO llal . at a gravitational wave frequency of 
/gw ~ HO (M/4OM ) _1 Hz. The dashed line shows the 
range of the interferometers for matched filtering performed 
with the EOB waveforms, which predict the evolution of the 
system up to the light ring orbit |0, Il9ll . at a gravitational 
wave frequency of /gw ~ 218 (M/40M©)" 1 Hz (notice that 
both these equations for /gw are for binaries of equal com- 
ponent masses). Since the EOB waveforms extend beyond 
the ISCO, they have longer duration and greater energy in the 
LIGO band which explains why the range for the EOB wave- 
forms is greater than the range for the TaylorT3 or Pade wave- 
forms (calculations performed with the Pade waveforms result 
in ranges similar to those given by the TaylorT3 waveforms). 
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FIG. 1: Range (distance at which an optimally oriented inspiraling 
binary of given total mass would produce a signal-to-noise ratio of 8) 
of the LIGO interferometers during S2. The error bars are calculated 
from the fluctuations of the noise in the LIGO interferometers during 
S2. 

During S2 the LI interferometer was the most sensitive 
with a range of 7 Mpc for a 10 — 10 A/ binary (calculated 
using the TaylorT3 waveform). However, since for the search 
described in this paper we demanded that our candidate events 
are seen in coincidence between the two LIGO observatories 
(as described in Sec. IV B> the overall range of the search was 
determined by the less sensitive LHO interferometers and thus 
was smaller than this maximum. 



IV. FILTERING 

A. Detection template family 

As was mentioned in Sec. [Ill] the gravitational wave sig- 
nal from inspiraling black hole binaries of high masses enters 
the LIGO frequency sensitivity band in the later stages, when 
the post-Newtonian approximation is beginning to lose valid- 
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ity and different versions of the approximation are beginning 
to substantially differ from each other. In order to detect these 
inspiral signals we need to use filters based on phenomeno- 
logical waveforms (instead of model-based waveforms) that 
cover the function space spanned by different versions of the 
late-inspiral post-Newtonian approximation. 

It must be emphasized at this point that black hole binaries 
with small component masses (corresponding to total mass 
up to 10 Mq) enter the S2 LIGO sensitivity band at an early 
enough stage of the inspiral that the signal can be adequately 
approximated by the stationary phase approximation to the 
standard post-Newtonian approximation. For those binaries 
it is not necessary to use phenomenological templates for the 
matched filtering; the standard post-Newtonian waveforms 
can be used as in the search for BNS inspirals. However, us- 
ing the phenomenological waveforms for those binaries does 
not limit the efficiency of the search fP^l . In this search, in or- 
der to treat all black hole binaries uniformly, we chose to use 
the BCV templates with parameters that span the component 
mass range from 3 to 20 M©. 

The phenomenological templates introduced in ifTotl match 
very well most physical waveform models that have been sug- 
gested in the literature for BBH coalescences. Even though 
they are not derived by calculations based on a specific phys- 
ical model they are inspired by the standard post-Newtonian 
inspiral waveforms. In the frequency domain, they are 



h(f) = A(f)e 
where the amplitude A(f) is 



"/>(/) 



/>0, 



A{f) = f- 7 / 6 (l-af' 3 )6{U-f) 



(2) 



(3) 



coefficients equal to and simplify the phase to 

i>(J) = ^o+27r/io + r 5/3 W>o+V'3/) (5) 
= 4o+Mf) (6) 

where the subscript s stands for "simplified". 

For the filtering of the data, a bank of BCV templates was 
constructed over the parameters / cut , i/jq and -03 (intrinsic tem- 
plate parameters). For details on how the templates in the 
bank were chosen see Sec. IV Bl For each template, the signal- 
to-noise ratio (defined in Sec. IIV B^ is maximized over the 
parameters to, 4>o and a (extrinsic template parameters). 



B. Filtering and signal-to-noise ratio maximization 

For a signal s, the signal-to-noise ratio (SNR) resulting 
from matched-filtering with a template h is 



p{h) 



(s,h) 



with the inner product (s, h) being 



(7) 



, h) _ , r s(f)h*(f) , f m r s(f)h*(f) 



and the phase tp{f) is 



and Sh(f) being the one-sided noise power spectral density. 

Various manipulations (given in detail in App. give the 
expression for the SNR (maximized over the extrinsic parame- 
ters (fio, a and to) that was used in this search. That expression 
is 



(4) 



n=0 



In Eq. ((3j 6> is the Heaviside step function and in Eq. (|4} to 

and (f>o are offsets on the time of arrival and on the phase of 
the signal respectively. Also, a, / cu t and i/j n are parameters 
of the phenomenological waveforms. 

Two components can be identified in the amplitude part 
of the BCV templates. The /~ 7 / 6 term comes from the 
restricted-Newtonian amplitude in the Stationary Phase Ap- 
proximation (SPA) lUllHEi. The term a/ 2 / 3 x J" 7 / 6 = 
af^ 1 / 2 is introduced to capture any post-Newtonian ampli- 
tude corrections and to give high overlaps between the BCV 
templates and the various models that evolve the binary past 
the ISCO frequency. Additionally, in order to obtain high 
matches with the various post-Newtonian models that predict 
different terminating frequencies, a cutoff frequency / cut is 
imposed to terminate the waveform. 

It has been shown 1 1 9fl that in order to achieve high matches 
with the various model-derived BBH inspiral waveforms it is 
in fact sufficient to use only the parameters tpo an d ^3 m the 
phase expression in Eq. 0, if those two parameters are al- 
lowed to take unphysical values. Thus, we set all other ip n 



^maximized - \ V 1*1 1 2 + 1^2 ? + ) (9) 

+ y\Fi\ 2 + \F2\ 2 -2Z(F 1 FZ). 



where 



df. (11) 



The quantities 01, 61 and 62 are dependent on the noise and the 
cutoff frequency / cut and are defined in App.|A| The original 
suggestion of Buonanno, Chen and Vallisneri was that for the 
SNR maximization over the parameter a the values of (a x 

fcit) should be restricted within the range [0, 1], for reasons 
that will be explained in Sec. lVCTI However, in order to be 
able to perform various investigations on the values of a we 
leave its value unconstrained in this maximization procedure. 
More details on this can be found in Sec. IV C Tl 
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V. SEARCH FOR EVENTS 
A. Data quality and veto study 

The matched filtering algorithm is optimal for data with 
a known calibrated noise spectrum that is Gaussian and 
stationary over the time scale of the data blocks analyzed 
(2048 s, described in Section IV B> . which requires stable, 
well-characterized interferometer performance. In practice, 
the performance is influenced by non- stationary optical align- 
ment, servo control settings, and environmental conditions. 
We used two strategies to avoid problematic data. The first 
strategy was to evaluate data quality over relatively long time 
intervals using several different tests. As in the BNS search, 
time intervals identified as being unsuitable for analysis were 
skipped when filtering the data. The second strategy was to 
look for signatures in environmental monitoring channels and 
auxiliary interferometer channels that would indicate an envi- 
ronmental disturbance or instrumental transient, allowing us 
to veto any candidate events recorded at that time. 

The most promising candidate for a veto channel was 
Ll:LSC-POB J (hereafter referred to as "POBI"), an auxiliary 
channel measuring signals proportional to the length fluctua- 
tions of the power recycling cavity. This channel was found 
to have highly variable noise at 70 Hz which coupled into the 
gravitational wave channel. Transients found in this channel 
were used as vetoes for the BNS search in the S2 data 0]. 
Hardware injections of simulated inspiral signals l23ll were 
used to prove that signals in POBI would not veto true inspi- 
ral gravitational waves present in the data. 

Investigations showed that using the correlations between 
POBI and the gravitational wave channel to veto candidate 
events would be less efficacious than it was in the BNS search. 
Therefore POBI was not used a an a-priori veto. However, the 
fact that correlations were proven to exist between the POBI 
signals and the BBH inspiral signals made it worthwhile to 
follow-up the BBH inspiral events that resulted at the end of 
our analysis and check if they were correlated with POBI sig- 
nals (see Sec. lVIIBl 

As in the BNS search in the S2 data, no instrumental vetoes 
were found for HI and H2. A more extensive discussion of 
the LIGO S2 binary inspiral veto studies can be found in |'24!i]. 



B. Analysis Pipeline 

In order to increase the confidence that a candidate event 
coming out of our analysis is a true gravitational wave and not 
due to environmental or instrumental noise we demanded the 
candidate event to be present in the LI interferometer and at 
least one of the LHO interferometers. Such an event would 
then be characterized as a potential inspiral event and be sub- 
ject to thorough examination. 

The analysis pipeline that was used to perform the BNS 
search (and was described in detail in |2|) was the starting 
structure for constructing the pipeline used in the BBH inspi- 
ral search described in this paper. However, due to the dif- 
ferent nature of the search, the details of some components of 



the pipeline needed to be modified. In order to highlight the 
differences of the two pipelines and to explain the reasons for 
those, we describe our pipeline below. 

First, various data quality cuts were applied on the data and 
the segments of good data for each interferometer were inden- 
tified. The times when each interferometer was in stable op- 
eration (called science segments) were used to construct three 
data sets corresponding to: (1) times when all three interfer- 
ometers were operating (L1-H1-H2 triple coincident data), (2) 
times when only the LI and HI (and not the H2) interfer- 
ometers were operating (LI -HI double coincident data) and 
(3) times when only the LI and H2 (and not the HI) interfer- 
ometers were operating (L1-H2 double coincident data). The 
analysis pipeline produced a list of coincident triggers (times 
and template parameters for which the SNR threshold was ex- 
ceeded and all cuts mentioned below were passed) for each of 
the three data sets. 

The science segments were analyzed in blocks of 2048 s us- 
ing the FINDCHIRP implementation 12511 of matched filtering 
for inspiral signals in the LIGO Algorithm library l2fill . The 
original version of FINDCHIRP had been coded for the BNS 
search and thus had to be modified to allow filtering of the 
data with the BCV templates described in Sec. lIVI 

The data for each 2048 s block was first down-sampled 
from 16384 Hz to 4096 Hz. It was subsequently high-pass 
filtered at 90 Hz in the time domain and a low frequency cut- 
off of 100 Hz was imposed in the frequency domain. The 
instrumental response for the block was calculated using the 
average value of the calibration (measured every minute) over 
the duration of the block. 

The breaking up of each segment for power spectrum es- 
timation and for matched-filtering was identical to the BNS 
search |2] and is briefly mentioned here so that the terminol- 
ogy is established for the pipeline description that follows. 

Triggers were not searched for within the first and last 64 s 
of a given block, so subsequent blocks were overlapped by 
128 s to ensure that all of the data in a continuous science seg- 
ment (except for the first and last 64 s) was searched. Any 
science segments shorter than 2048 s were ignored. If a sci- 
ence segment could not be exactly divided into overlapping 
blocks (as was usually the case) the remainder of the segment 
was covered by a special 2048 s block which overlapped with 
the previous block as much as necessary to allow it to reach 
the end of the segment. For this final block, a parameter was 
set to restrict the inspiral search to the time interval not cov- 
ered by any previous block, as shown in Fig.|2] 

Each block was further split into 15 analysis segments of 
length 256 s overlapped at the beginning and at the end by 64 
s. The average power spectrum Sh(f) for the 2048 s of data 
was estimated by taking the median of the power spectra of 
the 15 segments. We used the median instead of the mean to 
avoid biased estimates due to large outliers, produced by non- 
stationary data. The calibration was applied to the data in each 
analysis segment. 

In order to avoid end-effects due to wraparound of the dis- 
crete Fourier transform when performing the matched filter, 
the frequency-weighting factor 1/Sh(f) was truncated in the 
time domain so that its inverse Fourier transform had a max- 
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FIG. 2: The algorithm used to divide science segments into data anal- 
ysis segments. Science segments are divided into 2048 s blocks over- 
lapped by 128 s. (Science segments shorter than 2048 s are ignored.) 
An additional block with a larger overlap is added to cover any re- 
maining data at the end of a science segment. Each block is divided 
into 15 analysis sesgments of length 256 s for filtering. The first 
and last 64 s of each analysis segment are ignored, so the segments 
overlap by 128 s. Areas shaded black are searched for triggers by 
the search pipeline. The gray area in the last block of the science 
segment is not searched for triggers as this time is covered by the 
preceding block, although these data points are used in estimating 
the noise power spectral density for the final block. 



imum duration of ±16 s. The output of the matched filter 
near the beginning and end of each segment was corrupted 
by end-effects due to the finite duration of the power spec- 
trum weighting and the template. By ignoring the filter output 
within 64 s of the beginning and end of each segment, we 
ensured that only uncorrupted filter output is searched for in- 
spiral triggers. This necessitated the overlapping of segments 
and blocks described above. 

The single-sided power spectral density (PSD) of the noise 
Sh{f) in the LI interferometer was estimated independently 
for each LI block that was coincident with operation of at 
least one LHO interferometer. The PSD was used to construct 
a template bank for filtering that block. The bank was con- 
structed over the parameters xpo and xp3 so that there was no 
less than 95% overlap (defined in the sense of Eq. ( IA16i be- 
tween two neighboring (in {^0, ips, /cut} parameter space) 
templates, if the value of a of those templates was equal to 
0. The ipo — -03 space was tiled using a square grid based on 
the metric in Eq. (117) of fl9ll . For each pair of ipo an d ^3, 
three values of the cutoff frequency / cut were generated. This 
process gave the three intrinsic parameters for each template. 
Details on the exact values of the parameters used are given in 
Sec. IV C II The number of templates in the bank varied with 
the PSD. For this search the number of templates ranged be- 
tween 741 and 1296 templates per 2048 s LI analysis block, 
with the average number being 958. 

The data from the LI interferometer for the block was then 
matched-filtered, as described previously, against the bank of 
templates, with a SNR threshold p^ lcsh to produce a list of 
triggers. As will be explained below, the % 2 -veto ll27ll that 
was used for the BNS search was not used in this search. For 
each block in the LHO interferometers, a triggered bank was 



created consisting of every template that produced at least one 
trigger in the LI data during the time of the LHO block. This 
triggered bank was used to matched-filter the data from the 
LHO interferometers. For times when only the H2 interferom- 
eter was operating in coincidence with LI, the triggered bank 
was used to matched-filter the H2 blocks that overlapped with 
LI data. For all other times, all HI data that overlapped with 
LI data was matched-filtered using the triggered bank for that 
block. For HI triggers produced during times when all three 
interferometers were operating, a second triggered bank was 
produced for each H2 block. This triggered bank consisted of 
every template which produced at least one trigger found in 
coincidence in LI and HI during the time of the H2 block. 
The H2 block was matched-filtered with this bank. 

Before any triggers were tested for coincidence, all triggers 
with (a x / C 2 U { 3 ) greater than a threshold a* resh were rejected. 
The reason for this veto will be explained in Sec. IV CI 

For triggers to be considered coincident between two in- 
terferometers they had to be observed in both interferometers 
within a time window that allowed for the error in measure- 
ment of the time of the trigger. If the interferometers were not 
co-located, this parameter was increased by the light travel 
time between the two LIGO observatories (10 ms). We then 
ensured that the triggers had consistent waveform parameters 
by demanding that the two parameters ipo and ^3 for the tem- 
plate were exactly equal to each other. 

Triggers that were generated from the triple coincident data 
were required to be found in coincidence in the LI and HI 
interferometers. We searched the H2 data from these triple 
coincident times but did not reject Ll-Hl coincident triggers 
that were not found in the H2 data, even if, based on the SNR 
observed in HI they would be expected to be found in H2. 
This was a looser rejection algorithm than the one used for the 
BNS search | 2] and could potentially increase the number of 
false alarms in the triple coincident data. The reason for using 
this algorithm in the BBH search is explained in Sec. IV C 21 
where the coincidence parameter tuning is discussed. 

The last step of the pipeline was the clustering of the trig- 
gers. The clustering is necessary because both large astro- 
physical signals and instrumental noise bursts can produce 
many triggers with coalescence times within a few ms of each 
other and with different template parameters. Triggers sepa- 
rated by more than 0.25 s were considered distinct. This time 
was approximately half of the duration of the longest signal 
that we could detect in this search. We chose the trigger with 
the largest combined SNR from each cluster, where the com- 
bined SNR is defined in SeclVTl 

To perform the search on the full data set, a Directed 
Acyclic Graph (DAG) was constructed to describe the work 
flow, and execution of the pipeline tasks was managed by Con- 
dor 12811 on the various Beowulf clusters of the LIGO Scien- 
tific Collaboration. The software to perform all steps of the 
analysis pipeline and construct the DAG is available in the 
package LALAPPS lErjl . 
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C. Parameter Tuning 

An important part of the analysis was to decide on the val- 
ues of the various parameters of the search, such as the SNR 
thresholds and the coincidence parameters. The parameters 
were chosen so as to compromise between increasing the de- 
tection efficiency and lowering the number of false alarms. 

The tuning of all the parameters was done by studying the 
playground data only. In order to tune the parameters we per- 
formed a number of Monte-Carlo simulations, in which we 
added simulated BBH inspiral signals in the data and searched 
for them with our pipeline. While we used the phenomenolog- 
ical detection templates to perform the matched filtering, we 
used various model-based waveforms for the simulated sig- 
nals that we added in the data. Specifically, we chose to inject 
effective-one-body (EOB, 0,|1|1|1), Pade (PadeTL H) and 
standard post-Newtonian waveforms (TaylorT3, O), all of 
second post-Newtonian order. Injecting waveforms from dif- 
ferent families allowed us to additionally test the efficiency of 
the BCV templates for recovering signals predicted by differ- 
ent models. 

In contrast to neutron stars, there are no observation-based 
predictions about the population of BBH systems in the Uni- 
verse. For the purpose of tuning the parameters of our pipeline 
we decided to draw the signals to be added in the data from a 
population with distances between 10 kpc and 20 Mpc from 
the Earth. The random sky positions and orientations of the 
binaries resulted in some signals having much larger effec- 
tive distances (distance from which the binary would give the 
same signal in the data if it were optimally oriented). It was 
determined that using a uniform-distance or uniform-volume 
distribution for the binaries would overpopulate the larger dis- 
tances (for which the LIGO interferometers were not very sen- 
sitive during S2) and only give a small number of signals in 
the small-distance region, which would be insufficient for the 
parameter tuning. For that reason we decided to draw the sig- 
nals from a population that was uniform in log(distance). For 
the mass distribution, we limited each component mass be- 
tween 3 and 20 Mq. Populations with uniform distribution of 
total mass were injected for the tuning part of the analysis. 

There were two sets of parameters that we could tune in 
the pipeline: the single interferometer parameters, which were 
used in the matched filtering to generate inspiral triggers in 
each interferometer, and the parameters which were used to 
determine if triggers from different interferometers were co- 
incident. The single interferometer parameters that needed to 
be tuned were the ranges of values for ipo and ips in the tem- 
plate bank, the number of / cut frequencies for each pair of 
Wo, ^3} m the template bank and the SNR threshold p th ™ Bh . 
The coincidence parameters were the time coincidence win- 
dow for triggers from different interferometers, St, and the 
coincidence window for the template parameters ipo and -03. 
Due to the nature of the triggered search pipeline, parameter 
tuning was carried out in two stages. We first tuned the sin- 
gle interferometer parameters for the primary interferometer 
(LI). We then used the triggered template banks (generated 
from the L 1 triggers) to explore the single interferometer pa- 
rameters for the less sensitive LHO interferometers. Finally 



the parameters of the coincidence test were tuned. 

1. Single interferometer tuning 

Based on the playground injection analysis it was de- 
termined that the range of values for ipo had to be 
[10, 550000] Hz 5/3 and the range of values for ip 3 had to be 
[—4000, — 10] Hz 2 / 3 in order to have high detection efficiency 
for binaries of total mass between 6 and 40 Mq. 

Our numerical studies showed that using between 3 and 5 
cutoff frequencies per Wo, 1^3} pair would yield very high 
detection efficiency. Consideration of the computational cost 
of the search led us to use 3 cutoff frequencies per pair, 
thus reducing the number of templates in each bank by 40% 
compared to a template bank with 5 cutoff frequencies per 
Wo, iMpair. 

Our Monte-Carlo simulations showed that, in order to be 
able to distinguish an inspiral signal from an instrumental or 
environmental noise event in the data, the minimum require- 
ment should be that the trigger has SNR of at least 7 in each 
interferometer. A threshold of 6 (that was used in the BNS 
search) resulted in a very large number of noise triggers that 
needlessly complicated the data handling and post-pipeline 
processing. 

A standard part of the matched-filtering process is the x 2 - 
veto ll27ll . The x 2 - veto compares the SNR accumulated in 
each of a number of frequency bands of equal inspiral tem- 
plate power to the expected amount in each band. Gravita- 
tional waves from inspiraling binaries give small \ 2 values 
while instrumental artifacts give high \ 2 values. Thus, the 
triggers resulting from instrumental artifacts can be vetoed by 
requiring the value of x 2 for a trigger to be below a threshold. 
The test is very efficient at distinguishing BNS inspiral sig- 
nals from loud non-Gaussian noise events in the data and was 
used in the BNS inspiral search 01 in the S2 data. However, 
we found that the x 2 -veto was not suitable for the search for 
gravitational waves from BBH inspirals in the S2 data. The 
expected short duration, low bandwidth and small number of 
cycles in the S2 LIGO frequency band for many of the possi- 
ble BBH inspiral signals made such a test unreliable unless a 
very high threshold on the values of \ 2 were to be set. A high 
threshold, on the other hand, resulted in only a minimal re- 
duction in the number of noise events picked up. Additionally 
the x 2 -veto is computationally very costly. We thus decided 
to not use it in this search. 

As mentioned earlier, the SNR calculated using the BCV 
templates was maximized over the template parameter a. For 
every value of a, there is a frequency /o for which the ampli- 
tude factor (1 — a/ 2 / 3 ) becomes zero: 

fa = "~ 3/2 - (12) 

If the value of a associated with a trigger is such that the fre- 
quency /o is greater than the cutoff frequency / cut of the tem- 

2/3 

plate (and consequently a / c „ t < 1), then the high-frequency 
behavior of the phenomenological template is as expected for 
an inspiral gravitational waveform. If the value of a is such 
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FIG. 3: Time domain plots of BCV waveforms for different values 
of of- The top plot is for qf = 0, the middle is for oif = 1 and the 
bottom is for ct f = 2. For all three waveforms tpo = 150000 Hz 5 ' 3 , 
= -1500 Hz 2/3 and / cut = 500 Hz. It can be seen that the 
behavior is not that of a typical inspiral waveform for hf = 2. 

that /o is smaller than the cutoff frequency / cut (and conse- 
quently a / cl ( t > 1), the amplitude of the phenomenolog- 
ical waveform becomes zero before the cutoff frequency is 
reached. For such a waveform, the high-frequency behavior 
does not resemble that of a typical inspiral gravitational wave- 
form. For simplicity we define 

a F= a/cut 3 - (13) 

The behavior of the BCV waveforms for three different values 
of up is shown in Fig. [3] where it can be seen that for the case 
of ap = 2 the amplitude becomes zero and then increases 
again. 

Despite that fact, many of the simulated signals that we 
added in the data were in fact recovered with values of ap > 
1, with a higher SNR than the SNR they would have been re- 
covered with, had we imposed a restriction on a. Additionally 
some signals gave SNR smaller than the threshold for all val- 
ues of a that gave ap < 1. Multiple studies showed that this 
was due to the fact that we only had a limited number of cutoff 
frequencies in our template bank and in many cases the lack 
of the appropriate ending frequency was compensated for by 
a value of a that corresponded to an untypical inspiral gravi- 
tational waveform. 

We performed various investigations which showed that re- 
jecting triggers with ap > 1 allowed us to still have a very 
high efficiency in detecting BBH inspiral signals (although 
not as high as if we did not impose that cut) and the cut pri- 
marily affected signals that were recovered with SNR close to 
the threshold. It was also proven that such a cut reduced the 
number of noise triggers significantly, so that the false-alarm 
probability was significantly reduced as well. The result was 



that such a cut provided a clearer distinction between the noise 
triggers that resulted from our pipeline and the triggers that 
came from simulated signals injected in the data. In order to 
increase our confidence in the triggers that came out of the 
pipeline being BBH inspiral signals, we rejected all triggers 
with ap > 1 in this search. 

As mentioned in Sec. IIVI the initial suggestion of Buo- 
nanno, Chen and Vallisneri was that the parameter ap be con- 
strained from below to not take values less than zero. This 
suggestion was based on the fact that for values of a < 0, 

the amplitude factor (1 — afcut) can substantially deviate 
from the predictions of the post-Newtonian theory at high fre- 
quencies. Investigations similar to those described for the cut 
ap < 1 did not justify rejecting the triggers with ap < 0, so 
we set no low threshold for ap. 



2. Coincidence parameter tuning 

After the single interferometer parameters had been se- 
lected, the coincidence parameters were tuned using the trig- 
gers from the single interferometers. 

The time of arrival of a simulated signal at an interferome- 
ter could be measured within ±10 ms. Since the HI and H2 
interferometers are co-located, 5t was chosen to be 10 ms for 
H1-H2 coincidence. Since the light-travel time between the 
two LIGO observatories is 10 ms, St was chosen to be 20 ms 
for LHO-LLO coincidence. 

Because we performed a triggered search, the data from all 
three LIGO interferometers was filtered with the same tem- 
plates for each 2048 s segment. That led us to set the values 
for the template coincidence parameters Aipo and A ^3 equal 
to 0. We found that that was sufficient for the simulated BBH 
inspiral signals to be recovered in coincidence. 

The slight misalignment of the LI interferometer with re- 
spect to the LHO interferometers led us to choose to not im- 
pose an amplitude cut in triggers that came from the two dif- 
ferent observatories. This choice was identical to the choice 
made for the BNS search 01 ■ 

We considered imposing an amplitude cut on the triggers 
that came from triple coincident data and were otherwise co- 
incident between HI and H2. A similar cut was imposed on 
the equivalent triggers in the BNS search. The cut relied on 
the calculation of the "BNS range" for HI and H2. The BNS 
range is defined as the distance at which an optimally oriented 
neutron star binary, consisting of two components each of 1 .4 
Mq, would be detected with a SNR of 8 in the data. The value 
of the range depends on the PSD. For binary neutron stars the 
value of the range can be calculated for the 1 .4- 1 .4 Mq binary 
and then be rescaled for all masses. That is because the BNS 
inspiral signals always terminate at frequencies above 733 Hz 
(ISCO frequency for a 3-3 Mq binary, according to the test- 
mass approximation) and thus for BNS the larger part of the 
SNR comes from the high-sensitivity band of LIGO. For black 
hole binaries, on the other hand, the ending frequency of the 
inspiral varies from 1 10 Hz for a 20-20 Mq binary up to 733 
Hz for a 3-3 Mq binary (according to the test-mass approxi- 
mation). That means that the range depends not only on the 
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PSD but also on the binary that is used to calculate it. That 
can make a cut based on the range very unreliable and force 
rejection of triggers that should not be rejected. In order to 
be sure that we would not miss any BBH inspiral signals, we 
decided to not impose the amplitude consistency cut between 
HI and H2. 



VI. BACKGROUND ESTIMATION 

We estimated the rate of accidental coincidences (also 
known as background rate) for this search by introducing 
an artificial time shift At to the triggers coming from the 
LI interferometer relative to the LHO interferometers. The 
time-shift triggers were fed into the coincidence steps of the 
pipeline and, for the triple coincident data, to the step of the 
filtering of the H2 data and the H1-H2 coincidence. By choos- 
ing a shift larger than 20 ms (the time coincidence window 
between the two observatories), we ensured that a true gravi- 
tational wave could never produce coincident triggers in the 
time-shifted data streams. To avoid correlations, we used 
shifts longer than the duration of the longest waveform that 
we could detect (0.607 s given the low-frequency cutoff im- 
posed, as explained in Sec. lIIII . We chose to not time-shift the 
data from the two LHO interferometers relative to one another 
since there could be true correlations producing accidental co- 
incidence triggers due to environmental disturbances affecting 
both of them. The resulting time-shift triggers corresponded 
only to accidental coincidences of noise triggers. For a given 
time shift, the triggers that emerged from the pipeline were 
considered as one single trial representation of an output from 
a search if no signals were present in the data. 

A total of 80 time-shifts were performed and analyzed in 
order to estimate the background. The time shifts ranged from 
At = —407 s up to At = +407 s in increments of 10 s. The 
time shifts of ±7s were not performed. 



A. Distribution of background events 

Fig. |4] shows the time shift triggers that resulted from our 
pipeline (crosses) in the LHO SNR (pa) versus the LLO SNR 
(Pl) plane. There were only double coincident triggers in 
both the double and triple coincident data. Specifically, all 
the triggers present in the triple coincident data were Ll-Hl 
coincident and were not seen in H2. Thus, pn is defined as 
the SNR of either HI or H2, depending on which of the three 
S2 data sets described in Sec. IV Bl the trigger came from. 

There is a group of triggers at the lower left corner of the 
plot, which correspond to coincidences with SNR of no more 
than 15 in each observatory. There are also long "tails" of 
triggers which have high SNR (above 15) in one of the inter- 
ferometers and low SNR (below 10) in the other. The distribu- 
tion is quite different from the equivalent distribution that was 
observed in the BNS search 1 2|. The presence of the tails (and 
their absence from the corresponding distribution in the BNS 
search) can be attributed to the fact that the x 2 -veto was not 
applied in this search and thus some of the loud noise events 
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FIG. 4: The LLO and LHO SNRs of the accidental coincidences 
from the time-shifted triggers (crosses) and the triggers from the sim- 
ulated signal injections (circles) are shown. The dashed lines show 
the equal false-alarm contours. 



that would have been eliminated by that test have instead sur- 
vived. 

For comparison, the triggers from some recovered injected 
signals are also plotted in Fig. @] (circles). The distribution 
of those triggers is quite different from the distribution of 
the triggers resulting from accidental coincidences of noise 
events. The noise triggers are concentrated along the two axes 
of the (pl, pn) plane and the injection triggers are spread 
in the region below the equal-SNR diagonal line of the same 
plane. This distribution of the injection triggers is due to the 
fact that during the second science run the LI interferome- 
ter was more sensitive than the LHO interferometers. Con- 
sequently, a true gravitational wave signal that had compara- 
ble LLO and LHO effective distances would be observed with 
higher SNR in LLO than in LHO. The few injections that pro- 
duced triggers above the diagonal of the graph correspond to 
BBH systems that are better oriented for the LHO interferom- 
eters than for the LI interferometer, and thus have higher SNR 
at LHO than at LLO. 



B. Combined SNR 

We define a "combined SNR" for the coincident triggers 
that come out of the time-shift analysis. The combined SNR 
is a statistic based on the accidental coincidences and is de- 
fined so that the higher it is for a trigger, the less likely it is 
that the trigger is due to an accidental coincidence of single- 
interferometer uncorrected noise triggers. Looking at the plot 
of the pn versus the pl of the background triggers we notice 
that the appropriate contours for the triggers at the lower left 
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corner of the plot are concentric circles with the center at the 
origin. However, for the tails along the axes the appropriate 
contours are "L" shaped. The combination of those two kinds 
of contours gives the contours plotted with dashed lines in 
Fig-El Based on these contours, we define the combined SNR 
of a trigger to be 



pc = mm{W pi + p\, 2pH - 3, 2p L - 3}. 



(14) 



After the combined SNR is assigned to each pair of triggers, 
the triggers are clustered by keeping the one with the high- 
est combined SNR within 0.25 s, thus keeping the "highest 
confidence" trigger for each event. 



VII. RESULTS 

In this section we present the results of the search in the 
S2 data with the pipeline described in Sec. IV Bl The com- 
bined SNR was assigned to the candidate events according to 
Eq. (O. 



A. Comparison of the mishitted triggers to the background 

There were 25 distinct candidate events that survived all the 
analysis cuts. Of those, 7 were in the Ll-Hl double coincident 
data, 10 were in the L1-H2 double coincident data and 8 were 
in the L1-H1-H2 triple coincident data. Those 8 events ap- 
peared only in the LI and HI data streams and even though 
they were not seen in H2 they were still kept, according to the 
procedure described in Sec. IV Bl 
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FIG. 5: Expected accidental coincidences per S2 (triangles) with 
one standard deviation bars. The number of events in the S2 (circles) 
is overlayed. 

In order to determine if there was an excess of candidate 
events above the background in the S2 data, we compared the 



number of zero-shift events to the expected number of acci- 
dental coincidences in S2, as predicted by the time-shift anal- 
ysis described in Sec. IVII Fig. [5] shows the mean cumula- 
tive number of accidental coincidences (triangles) versus the 
combined SNR squared of those accidental coincidences. The 
bars indicate one standard deviation. The cumulative number 
of candidate events in the zero-shift S2 data is overlayed (cir- 
cles). It is clear that the candidate events are consistent with 
the background. 



B. Investigations of the zero-shift candidate events 

Even though the zero-shift candidate events are consistent 
with the background, we investigated them carefully. We first 
looked at the possibility of those candidate events being cor- 
related with events in the POBI channel, for the reasons de- 
scribed in Sec. IV Al 

It was determined that the loudest candidate event and three 
of the remaining candidate events that resulted from our anal- 
ysis were coincident with noise transients in POBI. That led 
us to believe that the source of these candidate events was in- 
strumental and that they were not due to gravitational waves. 
The rest of the candidate events were indistinguishable from 
the background events. 



C. Results of the Monte-Carlo simulations 

As mentioned in Sec. IV Bl the Monte-Carlo simulations al- 
low us not only to tune the parameters of the pipeline, but also 
to measure the efficiency of our search method. In this section 
we look in detail at the results of Monte-Carlo simulations in 
the full data set of the second science run. 

Due to the lack of observation-based predictions for the 
population of BBH systems in the Universe, the inspi- 
ral signals that we injected were uniformly distributed in 
log(distance), with distance varying from 10 kpc to 20 Mpc 
and uniformly distributed in component mass (this mass dis- 
tribution was proposed by 12911 ). with each component mass 
varying between 3 and 20 M Q . 

Fig- H] shows the efficiency of recovering the injected sig- 
nals (number of found injections of a given effective distance 
divided by the total number of injections of that effective dis- 
tance) versus the injected LHO-effective distance. We chose 
to plot the efficiency versus the LHO-effective distance rather 
than versus the LLO-effective distance since HI and H2 were 
less sensitive than LI during S2. 

Our analysis method had efficiency of at least 90% for re- 
covering BBH inspiral signals with LHO-effective distance 
less than 1 Mpc for the mass range we were exploring. It 
should be noted how the efficiency of our pipeline varied for 
different injected waveforms. It is clear from Fig. [6] that the 
efficiency for recovering EOB waveforms was higher than that 
for TaylorT3 or PadeTl waveforms for all distances. This 
is expected because the EOB waveforms have more power 
(longer duration and larger number of cycles) in the LIGO 
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FIG. 6: The efficiency versus the LHO-effective distance for the dif- 
ferent families of injected waveforms is shown. The dashed line rep- 
resents the efficiency for PadeTl injections, the dotted line represents 
the efficiency for standard post-Newtonian time-domain waveforms 
and the solid line represents the efficiency for effective-one-body 
waveforms. All injected waveforms were of second post-Newtonian 
order. Binomial error bars are shown. 



frequency band compared to the PadeTl and TaylorT3 wave- 
forms. 

Even though the main determining factor of whether a sim- 
ulated signal is recovered or not is its effective distance in the 
least sensitive interferometer, it is worth investigating the ef- 
ficiency of recovering injections as a function of the injected 
total mass. In order to limit the effect of the distance in the 
efficiency, we chose all the injections with LHO-effective dis- 
tance between 1 and 3 Mpc and plotted the efficiency versus 
total mass in Fig.0 The plot shows a decrease in the detection 
efficiency as the total mass increases, for the TaylorT3 and the 
PadeTl waveforms, but not for the EOB waveforms. The rea- 
son for this is the same as the one mentioned previously. For 
the TaylorT3 and the PadeTl waveforms, the higher the total 
mass of the binary system, the fewer the cycles in the LIGO 
band. That leads to the reduction in efficiency for those wave- 
forms. The EOB waveforms, on the other hand, extend be- 
yond the ISCO and thus have more cycles in the LIGO band 
than the TaylorT3 or the PadieTl waveforms. The power of 
the signal in band is sufficient for those waveforms to be de- 
tected with an approximately equal probability for the higher 
masses. 



VIII. UPPER LIMIT ON THE RATE OF BBH INSPIRALS 

As was mentioned previously, a reliable upper limit cannot 
be calculated for the rate of BBH inspirals as was calculated 
for the rate of BNS inspirals. The reason for this is two-fold. 
Firstly, the characteristics of the BBH population (such as spa- 
tial, mass and spin distributions) are not known, since no BBH 



FIG. 7: The efficiency versus the injected total mass of all injected 
signals with LHO-effective distance between 1 and 3 Mpc for the dif- 
ferent families of injected waveforms is shown. The dashed line rep- 
resents the efficiency for PadeTl injections, the dotted line represents 
the efficiency for standard post-Newtonian time-domain waveforms 
and the solid line represents the efficiency for effective-one-body 
waveforms. All injected waveforms were of second post-Newtonian 
order. Binomial error bars are shown. 



systems have ever been observed. In addition, the BCV tem- 
plates are not guaranteed to have a good overlap with the true 
BBH inspiral gravitational wave signals. However, work by 
various groups has given insights on the possible spatial dis- 
tribution and mass distribution of BBH systems 12911 . Assum- 
ing that the model-based inspiral waveforms proposed in the 
literature have good overlap with a true inspiral gravitational 
wave signal and because the BCV templates have a good over- 
lap with those model-based BBH inspiral waveforms, we used 
some of those predictions to give an interpretation of the result 
of our search. 



A. Upper limit calculation 

Following the notation used in 13011 . let 7Z indicate the 
rate of BBH coalescences per year per Milky Way Equiv- 
alent Galaxy (MWEG) and N G (p* ) indicate the number of 
MWEGs which our search probes at p > p*. The probability 
of observing an inspiral signal with p > p* in an observation 
time T is 

P{p> p*;K) = l-e- nTN ^P*\ (15) 

A trigger can arise from either an inspiral signal in the data or 
from background. If denotes the probability that all back- 
ground triggers have SNR less than p* , then the probability of 
observing one or more triggers with p > p* is 

P(p> p*;TZ,b) = l-P b e- nTNG ^">. (16) 
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Given the probability Pi,, the total observation time T, the 
SNR of the loudest event /9 max , and the number of MWEGs 
^Vc(Pmax) to which the search is sensitive, we find that a fre- 
quentist upper limit, at 90% confidence level, is 



n<n. 



2.303 + In P b 

TN G (p max ) 



(17) 



For 1Z > 1Zgo%, there is more than 90% probability that at 
least one event would be observed with SNR greater than 
Pmax- Details of this method of determining an upper limit 
can be found in 13111 . In particular, one obtains a conservative 
upper limit by setting Pb = 1. We adopt this approach below 
because of uncertainties in our background estimate. 

During the 350.4 hours of non-playground data used in 
this search, the highest combined SNR that was observed 
was 16.056. The number Nq can be calculated, as in the 
BNS search 01, using the Monte-Carlo simulations that were 
performed. The difference from the BNS search is that 
the injected signals were not drawn from an astrophysical 
population, but from a population that assumes a uniform 
log(distance) distribution. The way this difference is handled 
is explained below. 

Our model for the BBH population carried the following 
assumptions: 

(1) Black holes of mass between 3 and 20 M Q result exclu- 
sively from the evolution of stars, so that our BBH sources are 
present only in galaxies. 

(2) The field population of BBHs is distributed amongst 
galaxies in proportion to the galaxies' blue light (as was as- 
sumed for BNS systems 0]). 

(3) The component mass distribution is uniform, with val- 
ues ranging from 3 to 20 M Q I29H32II . 

(4) The component spins are negligibly small. 

(5) The waveforms are an equal mixture of EOB, PadeTl 
and TaylorT3 waveforms. 

(6) The sidereal times of the coalescences are distributed 
uniformly throughout the S2 run. 

Assumption (4) was made because the BCV templates used 
in this search were not intended to capture the amplitude 
modulations of the gravitational waveforms expected to re- 
sult from BBH systems with spinning components. However, 
studies performed by BCV 13311 showed that the BCV tem- 
plates do have high overlaps (90% on average, the average 
taken over one thousand initial spin orientations) with wave- 
forms of spinning BBH systems of comparable component 
masses. Templates that are more suitable for detection of spin- 
ning BBH have been developed l33ll (BCV2) and are being 
used in a search for the inspiral of such binaries in the S3 
LIGO data. 

Assumption (5) was probably the most ad-hoc assumption 
in our upper limit calculation. Since this calculation was pri- 
marily intended to be illustrative of how our results can be 
used to set an upper limit, the mix of the waveforms was cho- 
sen for simplicity. It should be apparent how to modify the 
calculation to fit a different population model. 

With assumptions (1) and (2) we determined the number 
of MWEGs in each logarithmic bin of LHO-effective distance 
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FIG. 8: The number of Milky Way equivalent galaxies (crosses) and 
the efficiency of the search (circles) as calculated by Eq. 1 1 91 ver- 
sus the LHO-effective distance (in Mpc) are shown. The cumulative 
number of MWEGs JVg (diamonds) versus LHO effective distance 
is overlayed. The horizontal axis has a logarithmic scale, in accor- 
dance with the uniform-log distance distribution of the injected BBH 
inspiral signals. No units are given for the vertical axis because it cor- 
responds to three different quantities plotted against LHO-effective 
distance. 



(Nq c , crosses in Fig. [8}. We used our Monte-Carlo simula- 
tions to determine the efficiency for detection of a source at 
LHO-effective distance d with combined SNR pc > Pmax, 
Pmax being the combined SNR of the loudest event observed 
in our search. Specifically, we define 



/(rf;Pmax) 
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Nf(d;p max ) 



Mnj(rf) 
Nf(d] /9 max ) 



Ninj(d) 



TaylorT3 
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^((ii/Pmax) 



N< 



inj(d) 



(18) 



PTl 



and that is also plotted in Fig. [8] (circles). The efficiencies 
for each waveform family individually are given in Table [Q 
Finally, AG(pmax) was calculated as 

oo 

iV G (Pmax) = f( d ' P*»*) X N G( d )- ( 19 ) 

Evaluating that sum we obtained Nq = 1.6603. Using 
Eq. O we obtained TZ 90 o /o = 35 year" 1 MWEG" 1 . 



B. Error analysis 

A detailed error analysis is necessary and was carefully 
done for the rate upper limit calculated for BNS inspirals in 
|2j]. For the upper limit calculated here, such a detailed analy- 
sis was not possible due to the lack of a reliable BBH popula- 
tion. However, we estimated the errors coming from calibra- 
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TABLE I: Efficiencies of recovering simulated BBH inspiral signals from three different waveform families for each LHO-effective distance 
bin. 
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tion uncertainties and the errors due to the limited number of 
injections performed. 

The effect of the calibration uncertainties was calculated as 
was done in Sec. IX A 2 of 01 . In principle, those uncertain- 
ties affect the combined SNR pc as 

Sp c <max{[4(^L) 2 +4(^H) 2 l 1/2 , 2(<t>h), 2(Sp L )\ 



Pc 



(20) 



where we modified Eq. (23) of 01 based on our Eq. dl4> for 
the combined SNR. However, for the calculation of the effect 
of this error on the rate upper limit, we were interested in how 
the calibration uncertainties would affect the combined SNR 
of the loudest event. Careful examination of the combined 
SNR of the loudest event showed that for that event the min- 
imum of the three possible values in Eq. dl4> was the value 
+ Ph) 1 ^ 2 ' so we calculated the error due to calibration 
uncertainties by 



Sp c < 



4(^ L ) 2 

Pc 



Pc 



1/2 



(21) 



We simplified the calculation by being more conservative and 
using 



Spc< (Spl) 2 + (Sp H f 



1/2 



(22) 



We also used the fact that the maximum calibration errors at 
each site were 8.5% for LI and 4.5% for LHO (as explained 



in 1 21) to obtain 

Spc < [(0.085 pr x ) 2 + (0.045p£ ax ) 2 
The resulting error in N G was 

SN G \ ca i = ±0.0859 MWEG. 



1/2 



(23) 



(24) 



The errors due to the limited number of injections in our 
Monte-Carlo simulations had to be calculated for each loga- 
rithmic distance bin and the resulting errors to be combined in 
quadrature. Specifically, 



SN G \ 



MC 



oo 

{E[w 



x N G c (d) 



2-1 V 2 



(25) 



Because / was calculated using the 3 different waveform fam- 
ilies, the error Sf(d; p max ) is 

v rx- N f (d; Pmax )(N ini (d) - A^/w))! 1 / 2 

Sj{d;p„ 



3 2 [N ini (dW 



(26) 

where the sum was calculated over the three waveform fami- 
lies: EOB, PadeTl and TaylorT3. This gave 

SN G \ MC = ±0.0211 MWEG. (27) 

Finally we added the errors in Ag in quadrature and ob- 
tained 



SNr 



±0.0885 MWEG. 



(28) 
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Both contributions to this error can be thought of a 1-er vari- 
ations. In order to calculate the 90% level of the systematic 
errors we multiplied SNq by 1 .6, so 



SN, 



G\90% 



= ±0.1415 MWEG. 



(29) 



To be conservative, we assumed a downward excursion Nq = 
(1.6603 - 0.1415) MWEG = 1.5188 MWEG. When substi- 
tuted into the rate upper limit equation this gave 

TZ 90% = 38 year" 1 MWEG" 1 . (30) 



IX. CONCLUSIONS AND FUTURE PROSPECTS 

We performed the first search for binary black hole inspiral 
signals in data from the LIGO interferometers. This search, 
even though similar in some ways to the binary neutron star 
inspiral search, has some significant differences and presents 
unique challenges. There were no events that could be inden- 
tified as gravitational waves. 

The fact that the performance and sensitivity of the LIGO 
interferometers is improving and the frequency sensitivity 
band is being extended to lower frequencies makes us hopeful 
that the first detection of gravitational waves from the inspi- 
ral phase of binary black hole coalescences may happen in 
the near future. In the absense of a detection, astrophysically 
interesting results can be expected by LIGO very soon. The 
current most optimistic rates for BBH coalescences are of the 



order of 10" 4 year -1 MWEG" yj. It is estimated that at 
design sensitivity the LIGO interferometers will be able to de- 
tect binary black hole inspirals in at least 5600 MWEGs with 
the most optimistic calculations giving up to 13600 MWEGs 
l35ll . A science run of two years at design sensitivity is ex- 
pected to give BBH coalescence rate upper limits of less than 
10" 4 year" 1 MWEG -1 . 



APPENDIX A: FILTERING DETAILS 

The amplitude part of the BCV templates h(f) can be de- 
composed into two pieces, which are linear combinations of 
/ -7 / 6 and J -1 / 2 . Those expressions can be used to construct 
an orthogonal basis hj for the 4-dimensional linear subspace 
of templates with 4>q £ [0, 2tt) and a G (— oo, +oo). Specifi- 
cally, we want the basis vectors to satisfy 



h k , hj 



(Al) 



To do this we construct two real functions Ai(f) and A 2 (f), 
linear combinations of /~ 7 / 6 and / -1 / 2 , which are related to 
the four basis vectors via: 



k, 2 (f) = MM)^ U) 

hAf) = A lft {f)i^' U) - 
Then, Eq. IAH becomes: 

Mf)Mf). 



s h (j) 
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(A2) 
(A3) 



(A4) 



Since the templates have to be normalized, A\(f) and ^4.2 (/) 
must satisfy 
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Imposing condition JA4t gives the numerical values of the 
normalization factors. Those are 



«i 



-1/2 



7/3 



T 2 

1 5/3\-l/ 2 



'7/3 



where the integrals I k are 



4=4 



1 7/3 
J 2 

-'5/3 N-1/2 



df 



o f k S h (f)' 



(A5) 
(A6) 

(A7) 
(A8) 

(A9) 



The next step is to write the normalized template in terms 
of the 4 basis vectors 

Hf) = cihi(f) + c 2 h 2 (f) + c 3 h 3 (f) + cMl) (A10) 
with 



Cl 

c 2 

C;', 
C4 

where uj is related to a by 

tanoj : 



cos 0o cosw, 
cos </>o sin w, 
sin </>o cosw, 
sin <fio sinw, 



b 2 + ha 



(All) 
(A12) 
(A13) 
(A14) 



(A15) 



Once the filters are designed, the overlap is calculated and 
is equal to 



p= {s,h) = K\ cos ui cos 0o + K 2 sin uj cos </>o(A16) 
+K 3 cos w sin </> + K4 sin u> sin O 



where Kj = (s,hj),k= 1, 2, 3, 4 are the four integrals that 
are necessary, namely 

45(/) ai /-V« „_^ fn , 



A'i = 3? 



S h (f) 



'df, 



(A17) 



A" 2 = 3? 



/cut 45(/)(6 1 /- 7 / 6 + b 2 f- 1/2 ) c ^ s(f)dL 



S h (f) 



(A18) 



k, = _ a tSMflsW* (A19) 
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Ea. lA22l can then be written as 



Ka 



fc 



4 ±s(f)(hr 7/6 + b 2 f- 1/2 ) 



-iips(f) 



/O -W) 

Maximizing the SNR over <fro and oj we get 



(A20) 



P 



maximize 



d 



y^i+A' 4 ) 2 + (^-A' 3 ) 2 (A21) 



+ i^ 1 -A^ + (if 2 + if 3 ) 2 . 
The values of 0o and a that give the maximized SNR are 



1 



. arctan ^±^-i arctan (A22) 
K X -K A 2 Ai+if 4 



-'maximize' 



zed 



\F 2 1 2 + 23(FiF 2 *) (A26) 



29(FiF 2 *). 
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&2 tan w E 



ai + b\ tan w max ' 

where 

^ = \ arctan + i arctan . (A24) 

Z Ai + A4 Z Ai — A4 

An extensive discussion on the values of a is provided in 
Sec. lVCTl 

An equivalent expression for the SNR, which is computa- 
tionally less costly, can be produced if we define 



F 1 = Ki - iK 3 , 
F 2 = K 2 - %K A . 
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